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Abstract. We study the spherical mean transform on I". The transform 
is characterized by the Euler-Poisson-Darboux equation. By looking at the 
spherical harmonic expansions, we obtain a system of 1 + 1-dimension hy- 
perbolic equations, which provide a good machinery to attack problems of 
spherical mean transform. 

As showcases, we discuss two known problems. The first one is a local 
uniqueness problem investigated by M. Agranovsky and P. Kuchment, [Mem- 
oirs on Differential Equations and Mathematical Physics, 52:1—16, 2011]. We 
present a simple proof which works even under a weaker condition. The second 
problem is to characterize the kernel of spherical mean transform on annular 
regions, which was studied by C. Epstein and B. Kleiner [Comm. Pure Appl. 
Math., 46(3):441— 451, 1993]. We present a short proof that simultaneously ob- 
tains the necessity and sufficiency for the characterization. As a consequence, 
we derive a reconstruction procedure for the transform with additional interior 
(or exterior) information. 

We also discuss how the approach works for the hyperbolic and spherical 
spaces. 



Let / be a function defined on R n . The spherical mean transform 11(f) of / is 
defined as 



Here, S(x,t) C K™ is the sphere of radius t centered at x, da(x) is the surface 
measure on S(x,t), and \S(x,t)\ is the total measure of S(x,t). The spherical 
mean transform has been intensively studied due to its applications to PDEs, 
approximation theory, inverse scattering, and thermoacoustic tomography (e.g., 



[CH621 iJohSll |Bey83a[ |Biy83bl IFHR071 ILTM1 ILTMl |AQ96l ) . 



The spherical mean transform can be characterized by the Darboux (or Euler- 
Poisson-Darboux) equation. Namely, if / is smooth enough (for example / 6 



C°°(R n )), then G(x,t) = K(f)(x,t) satisfies (e.g., (Joh81] 1: 

m f G tt (x,t) + B=iG t (x,t) - A x G(x,t) = 0, (x,t) e R™ x R + , 

[ > \ G(x,0) = f(x), G t (x,0) =0, xeR n . 

Conversely, if G(x,t) 6 C°°(IR™ x R + ) satisfies the above equation then G = 11(f). 
Let us recall the polar coordinate decomposition of the Laplacian in R": 
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S(x,t) 
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n — 1 d 
r dr 



+ — A s „-i, 
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where r = \x\ and A§n-i is the Laplacc-Beltrami operator on the unit sphere S n 1 . 
Hence, for any twice diffcrentiable function g = g(r): 



A[g(r)r m Y m {8)} = r m Y m (6) 



d n — 1 + 2m d 
dr 2 r dr 



9(r), 



where Y m is any spherical harmonics of degree m. 

Let us expand / and G in terms of spherical harmonics: 



oo /„ 



( 2 ) /(*) = EE^W r ^ m (^ G(x,t) = Y,Y.9ml{r,t)r m Yr{e). 



m=0 1 = 1 k=0 1 = 1 



For each m > 0, we denote B r ,-, 
(3) 



d , n-l+2m 
dr' 2 r i 



Then, ([J} is equivalent to 



Bt,og m i(r,t) - B r , m g ml {r,t) = 0, (r,t) e K x R + , 
3mi(r,0) = / m ;(r), d t g m i(r,0) = 0, r e R, 

for each m > and I = 1, .., 2 m . This equation can be symmetrized by a differential 
operator. Indeed, let 

Q — ( T ^ + l) ( r ® + 1\ ( r ® + 

m \n dr J \n + 2dr J \n + 2(m ~ 1) dr 
We then have (see Appendix): 

(4) QmBr^n = B r ,oQm- 

Applying Q m to ([3]), we obtain the symmetric equation 

(5) [B t ,o-iB rj o]a(r,t)=0, 

where a = Q m gmi, for any m G N and I = 1, .., l m . 

In this article, we exploit equations ([3]) and ([5]) to present some short and clear 
proofs for some (old and new) properties of spherical mean transform. 

The article is organized as follows. In Section we revisit a local uniqueness 
result by Agranovsky and Kuchment |AK11) . We also point out that the same 
result still holds under a weaker condition. We then discuss a result by Epstein and 
Kleiner |EK93[ for spherical mean transform on the annular region in Section [3] We 
provide a reconstruction procedure for the transform when some additional interior 
(or exterior) information is provided. Finally, we discuss how the decompostion 
and symmetrization work on the hyperbolic and spherical spaces. 



2. A LOCAL UNIQUENESS RESULT 

Let us recall the following result by Agranovsky and Kuchment [AKlli Theorem 
10]: 

Theorem 2.1. Let f G C°° {B r +e ) such that f(x) = for x £ Br. Assume that 
lZ(f)(x, R) — for all x G B £ . Then, f(x) — for x G Br +£ . Here, Br is the 
sphere of radius R centered at 0. 

The above theorem was proved by utilizing the notion of ridge function and 
Titchmarsh theorem. We present here an alternative proof using equation ([3|) and 
simple energy arguments. Let us start with a basic domain of dependence argument 
which will be repeatedly used in this article: 
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Proposition 2.2. Let c > 0, e > 0, and u satisfy the equation: 

| + A v) »(r,V) = 0' V e K", r > e, 

1 u(e,y) =u e (y), u r (e,y) = 0, y G K™. 



(6) 



Assume that u e (y) = /or aZZ y 6 -B(j/Oi r o)- Then, u(r,y) = if r > e and 
\y-Vo\ +r < r +e. 

The proof Proposition 12.21 is well known (see, e.g., [CH62 ). It consists in only a 
simple energy argument. We will provide in Appendix for the sake of completeness. 

Proof of Theorem \2.1\ Let us consider the series expansion @ of / and G = TZ(f). 
We observe that the conditions f(x) — for x G Br and G(x, R) = lZ(f )(x, R) = 
for x G B £ imply 



(J) 



^(0,i) = 0, tG[0,iq, and g m i(r,R) = 0, s6[0,e] 



The goal is to prove that f m i(t) = g m i(0,t) = for t G [0, R + e}. Our idea is to 
transform equation @: 



i)! 1 



n-1 a 
t dt 



g m (r, t) - 



n— l+2ro 



5m (r,i) = 0, (r,t) G 



dr 2 r dr 

{ g m i (r, 0) = f rnl (r) , d t9ml (r, 0) = 0, Vr G M, 

to a Darboux equation for which r is the temporal variable. Indeed, let us introduce 
the function u(r,y) = gmiir, \y\), which is radially symmetric with respect to the 
variable y G K". Then for any y such that \y\ = t, 

d 2 n — 1 d 
dt 2 + t di_ 
We, hence, obtain the Darboux equation: 
d 2 n — 1 + 2m d 



9m(r,t). 



(8) 



u(r, y) - A y u(r, y) = 0. 



dr 2 r dr 

From (0, we obtain the following zero initial and boundary values of u: 

u(0, y) = 0, Vy G B R , u{r, y) = 0, V(r, y) G [0, e] X 5^. 

A simple energy argument shows that u(r, y) = for all (r, y) G i?« x [0, e]. Indeed, 
let 

£(r)= J [Mr,i/)| 2 + |Vu(r,y)| 2 ]dy. 

B(0,R) 

Since u(r, y) = for y G Sr, integration by parts gives 
dE(r) n 



dr 



[u r {r, y) - Au(r, y)] u r {r, y)dy. 



Due to ©, 



B(0,R) 
dE{r) 2(n - 1) 



dr 



\u r {r iy )\ 2 dy <0. 



B(0,iJ) 

We obtain E{r) < E(0) = for all < r < e. This implies E(r) = or u(r, y) = 
for all (r,y) e[0,e]x% 
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We now arrive to the equation 

f (& + *= 3 ^&-A w )«(r, !/)=<), 
\ u(s,y) = u r (e,y) = 0, y e M™. 

Due to Proposition 12. 2[ we obtain u(y,t) = for j/| + r < i? + e. This means 
gml( r ; t) = for all (r, t) such that r+£ < R+e. In particular, f m i(s) = g m i(0, s) = 
for s G [0, i? + e]. 

We conclude, due to the expansion ©, f(x) =0 for all x G Br +£ . □ 

The condition /(a;) = for i G -B_r is equivalent to TZ(f)(x,t) = for all (a?, t) 
such that \x\+t < R. It is much stronger than: D%H(f)(0, t) = for any multiindex 
a and i G [0, R}. Following the above proof, we obtain the same conclusion under 
this weaker condition: 

Theorem 2.3. Let f G C°°(B R+e ) such that D%U(f)(0, t) = for all multi-indices 
a and t G [0, R\. Assume that R(f){x, R) = for all x G B £ . Then fix) = for 
all x G -Br+ e . 

3. Spherical mean transform on the annular region 

Given < a < A, let us consider the annular region 

Ann(a, A) = {x e M" : a < \x\ < A}. 

We denote by Z°°(a,A) the space of all functions / G C°°(Ann(a, A)) satisfying 
Tl{f)(x,t) = for all (x,t) G A. Here, 

(9) A = {{x, t) : S(x, t) C Ann(a, A) and 5(0, a) C B(x, t)}, 

where B(x, t) is the ball of radius t centered at x. The following result characterize 
all the functions / G Z°°(a, A): 

Theorem 3.1. A function f belongs to Z°°{a 1 A) if and only if fo = and f m i is 
of the form 

k-l 

(10) /mK^E^™"^' 

i=0 

for any m > and I — 1, .., l m . 

This result was proved in [EK93|. The necessity was obtained by a projection 
formula on spaces of spherical harmonics. The sufficiency was proved by a series 
of lemmas concerning the connection between the functions in Z co {a, oo) and har- 
monic functions with certain behavior at infinity. In this section, we present an 
argument to simultaneously prove the necessity and sufficiency. It also provides a 
reconstruction procedure for the transform with additional interior/exterior infor- 
mation, which might find applications in biomedical imaging. 

Remark 3.2. By some convolution arguments as in EK93 , we can deduce from 
Theorem \3.1\ the same characterization for functions f G Z{a,A). Here, Z{a,A) is 
the set of functions f G C(Ann(a, A)) such that lZ(f)(x,t) = 0, \/(x,t) G A. 

Without loss of generality, we might also assume that / G C co (Ann(a, A)). 
Otherwise, we prove the same characterization on the region Ann{a + e, A — e) (for 
small enough e > 0) and then let e — > 0. 
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Let us extend / smoothly to W 1 . Then, G = 11(f) £ C°°(M™ x ' 
equation ([T]). The condition / £ Z oc (a,A) is equivalent to 



satisfies the 



(11) G(x,t) = 0, (x,t) e 1" x (a, A), such that |x| + 



a + A 



< 



A - a 



Lemma 3.3. Consider the decomposition of G as in fl|). Then f £ Z°°(a,A) if 
and only if 

(12) <?mi(0,*) = 0, We (a, A), 

for all m £ N and I = 1, .., l m . 

Proof. It suffices to prove that (fTTj) and (|T2|) are equivalent. The fact that (fTTj) 
implies (|12[) is obvious. We now prove the other implication. Similar to the proof 
of Theorem 12.11 we translate equation §3§ to a Darboux equation by introducing 
the function u(r,y) = g m i{'r 1 \y\). We then obtain: 

dr 2 



n — 1 + 2m d 
r dr 



u(r,y) = 0. 



Let 2/0 £ such that |y | = 2 -j^- The condition g m i(0, t) = for £ £ (a, A) implies 

"(0, 2/) = 0, for all y £ B(y , r ), 

where ro = . The domain of dependence argument (see Proposition 12.21) then 
implies 

u(r, y) = 0, for all (y, r) £ K" x R + such that |y — t/o| + T < tq. 
Recalling that u(r,y) = g m i{f, \y\), we obtain 

gmi(r,t) — 0, for all (x,t) such that \t — |?/o|| + t < tq. 

That is, 



gmi(r,t) = 0, for all (r, t) such that 



< - 



a - 



A 



+ r < 



A - a 



Since this is true for all m £ N and I = 1, .., l m , we conclude 

a + A 



G(x,t) = 0, for all (x,t) such that 
This finishes our proof. 



t - 



+ |a;| < 



A 



□ 



We now present our proof for Theorem 13.11 It simultaneously provides the 
necessity and sufficiency of the characterization in Theorem 13.11 



Proof of Theorem \'3.1\ Let a(r, t) = Q m gmi{ r i't), where 



(13) 



Qm = 



r d 
n dr 



i) 



2dr 







n + 2(m — 1) dr 



Due to Q , we arrive to the symmetric equation ([5]) : 



dt 2 



n — 1 d 



a(r,t) 



t dt 

We obtain (e.g., jJoh81l IHel84] 1 



r d 2 n — ld 
dr 2 



dr 



a{r,t)=0, (r,t) e 



a(s, 0) = a(0, s). 
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Simple observations show: 

a(0, s) =fl m j(0,s), a(0,s) = [Q m f m i](s). 

Therefore, 

(14) [QmfmlKs) =9ml{0,s). 

From Lemma [3T3l / G Z°°(a,A) if and only if gw(0, s) = for s G (a, A). Or, 
equivalently, 

(15) [em/ml](*)=0, VsG(o,A). 
Due to formula (fT5|) of Q m , (|15l) is equivalent to 



fml W 



E c^r""- 2 S ifro>l, 
0, if to = 0. 



This finishes our proof. □ 

We are now interested in reconstructing / on Ann(a, A) from 7t(f)\A- The prob- 
lem does not have a unique answer since the kernel described in Theorem 13.11 is 
nontrivial. We, therefore, make one more assumption: / G C°°(Ann(ro, A)) and 
it is known in the interior region Ann(ro, a] = {x G K™ : ro < \x\ < a}, for some 
number < ro < a. The problem of reconstructing / from TZ(f)\_A and f\Ann(ra,a] 
resembles interior tomography with prior interior knowledge, which is widely in- 
vestigated in biomedical imaging (e.g.. |KCND08| ICNDK08] ). The following result 
provide a formula useful for the reconstruction: 

Theorem 3.4. Consider the spherical harmonic decomposition of f and G = 
K(f). Let 

k ml (r) = " ( " + 2( 2 J:^! ( ^," 1)] r-'" +2 ^-^ /^Mr-V -r^dr. 

a 

Then for r G M.+ : 



m—1 

E 4nf- n ~ 2 \ ifm>l, 



(16) q m i{r) ■■= f m i(r) - k ml (r) = ^ 

0, if to = 
Proof. We claim that 

(17) Qmkml = 9 ml- 

Indeed, direct calculations show: 
r 



m-2)] 



n + 2(m — 



_J^ + 1 ) k ,( r ) _ n...[n + 2(m-2)} {n+2( 

r 

x [ g^O^y 1 - 1 ^ 2 -T 2 ) m - 2 dr 



We recall that 



, r d \ ( r d \ ( r d 

k fiW / \n + 2or / \rt + 2(m — 1) or 
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By induction, we obtain 



t 3 

Qrnklm(r) = I - •5- 



n dr 



g m i{0,T)r n - L dT 



9mi(0,r) 



This proves (fT"7|). Since / m ; also satisfies the same equation (|14[) . 

Qm(fml — 9ml) = 0. 

Therefore, g m /(r) := / m /(r) - k m i(r) has the form (JTHJ) . □ 
We now arrive to a procedure to compute / on Ann(a, A) from TZf\^ and 

•/*lArm(ro,a] ■ 

1) Compute f m i(r) for r G (r ,a] from /Unn(r 0l a]- 

2) Compute g m i(0, s) for s G {a, A) from G = lZ(f)\j\. 

3) Compute q m i{s) = fmi{s) — k m i{s) and its derivatives at s — a from the 
above knowledge of f m i and g m i . Use them to compute c l ml . 

4) Compute f m i = k m i + q m i on (a, A). Then, compute / on Ann(a, A). 

Remark 3.5. The same procedure also works when f(x) is provided in an exterior 
domain Ann[A, R) = {A < \x\ < R}, for some R > A. 

4. The approach for Hyperbolic and Spherical spaces 

We have considered the spherical mean transform from a PDE point of view. 
We decomposed the Euler-Poisson-Darboux equation into the 1 + 1-dimension hy- 
perbolic equations. By exploiting them and their symmetrized versions, we obtain 
some old and new properties for spherical mean transform. In this section, we 
describe the approach for the hyperbolic and spherical spaces. 

Let us consider the spherical mean transform on the hyperbolic space H n , which 
is the unit ball with the metric 

ds 2 = dx 2 

(l-|x| 2 ) 2 

The Laplace-Bcltrami operator A on H n is A = |(1 - |x| 2 )"V ((1 - |x| 2 ) 2 -"V) . 
In terms of polar coordinates, 

d 2 5 1 

(18) A = — + (n - 1) coth(r)— + A s »-i . 

dr 2 dr smlr (r) 

Here, r = du™(x,0) and A§ n -i is the Laplace-Beltrami operator on the Euclidean 
unit sphere S n_1 G R™ applying to the variable = t|t. 

The spherical mean transform IZ(f) of a function / is characterized by (e.g., 
[Hel84] ): 



(19) 



[&r + (n - 1) coth(t)f - A]G(x, t) = 0, (x, t) G 
G{x,0) = f{x), G t (x,0) = 0, xgHP. 



Conversely, if G(x,t) G C°°(lHI n x R + ) satisfies the above equation, then G(x,t) 
lZ(f)(x,t). Let us consider the decompositions^]: 

(20) f{x) = EE^M^W. G{x,t)=Y J Y,9mi{r,t)Y l m {0)- 

m=0 1=1 k=Q 1=1 



^This is different from the decomposition J2} on the Euclidean space. 
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Let us consider the operator 



d . d m(m + n — 2) 



V m ,r = ^-o + (" - 1) coth(r) 



2 



r 



Since Agn-iYJ" 1 — ~ni(m + n — 2)Y l m , for any function g = g(r): 

A[g(r)Yr(0)} = (V k g)(r)Yr(9). 
We obtain the following analog of equation ([3]) in Section [TJ 



(21) 



(X>0,t - T>m,r)9ml(r,t) = 0, 

flw(r,Q) = /wfr), §:g m i(r,0) = 0. 

This equation can be symmetrized by the operator Q m = ri...r m , where = 
+ (n + k — 2) coth(r). Indeed, one can verify (see |Ngull| ): 

Due to (|21[) . we obtain the symmetric equation for a(r,t) = Q m 3mi( r ;^) : 
(22) (V . t -V . r )a{r,t) =0. 

This equation provides a good tool to work with spherical mean transform on H™. 
For example, it was exploited in |Ngull| to characterize the functions / such that 
d£ft(/)(0, R) = for a fixed R > and all multi-index a. 

For the spherical space S™, which is the unit sphere in R™ +1 , the same argu- 
ment as above works; except that one need to replace the hyperbolic trigonometric 
functions (cosh, sinh, coth, etc) by the usual ones (cos, sin, cot, etc). 

Finally, we remark that one can use the same approach to investigate the ball 
transform (e.g., |Vol03] ) . instead of the spherical mean transform. 

Appendix 

In this Appendix, we present the proofs of equation |4] and Proposition 12.21 
Proof of equation [4j Let us recall that 

r a \ i r a \ I r o \ I r a 



i=0 



q TT ( + 1 ) ( + 1 ) ( + 1 

m J- J. \ r> _]_ 9i fir I \ n fir I \ n -1- 9 fir 



n + 2idr ) \ndr J \n + 2dr J \n + 2(m—l)dr 



One could easily check: 



d \ ( r d 

1 )B r , k = B r , k ^ — — — -- + 1 



n + 2(k-l)dr J ' ' \n + 2(k-l)dr 
By induction, we obtain: 

QmB r.m — •m • 

Proof of Proposition [2721 For e < r < tq + e, let 

E(r)= J [\u r (r,y)\ 2 + \Vu(r,y)\ 2 ]dy. 

B(y ,r +e—r) 



□ 
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Then 

dE(r) 



= 2 / [u r {r,y)u rr (r, y) + Vu(r,y)Vu r (r, y)]dy 
{\u r {r,y)\ 2 + \Vu(r,y)\ 2 )dy. 



dr 

B(yo,ro+e—r) 



1 



S(yo,ro+e-r) 

Taking integration by parts, we obtain 
dE(r 



= - I [u r (r,y)u rr (r,y) - Au{r,y)u r (r,y)]dy 

) 

[K(r, y)\ 2 - 2u r (r,y)d v u(r,y) + \Vu(r,y)\ 2 ] dy. 



dr 

B(yo,ro+e—r) 



S(yo,ra+e—r) 

Here, d v is the outer normal derivative on S(yo, r$ + e — r). Since \d u u\ < ||Vit|| 
the last integral is nonnegative. Therefore, 

dE(r) 



dr 



<2 / [u r {r,y)u rr (y,r) - Au(r,y)u r (r,y)]dy 

B(ya,r +6-r) 

n + 2m — 1 
rr(r, y) H u r [r, y) - Au(r, y) 



< 2 

B(Vo,ro+e— r) 



u r (r,y)dy. 



Due to equation (j6|, we obtain dE }p < for all e < r < r + e. Since E(e) = 0, 
we obtain E(r) = for all e < r < ro + e. This implies u(r, y) = for all 
y G B(yo,ro + e — r). Equivalently, u(r,y) — for all (r, y) such that r > e and 
\y ~ 2/o ] + r < r + e. □ 
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